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VORTEX  SHEDDING  AND  LASER  ABLATION 


The  Fayleigh-Taylor1  (RT)  instability  is  a  potential  obstacle  to  iner- 
tial  confinement  fusion  in  that  it  can  destroy  the  symmetry  of  the  inploding 
shell,  the  RT  instability  can  occur  anywhere  in  the  shell  when  the  profiles 
become  baroclinic,  i.e. ,  possess,  no  n-^ol  linear  density  and  pressure 
gradients,  the  ablation  layer,  which  separates  the  cold,  high-density  shell 
from  the  hot,  low-density  ablating  plasma,  is  particularly  susceptible  to  RT 
growth.  Several  researchers2' 3,  including  ourselves4' 5,  have  reported 
the  results  of  numerical  sinultions  of  the  RT  instability  in  laser  ablatively 
accelerated  targets,  the  results  were  consistent  only  in  that  the  linear 
growth  rates  were  below  the  classically  predicted  values  (y  *  (kg)l/2). 
Several  mechanisms  have  been  proposed  to  explain  the  reduced  linear  growth  of 
the  instability  including  thin  layer  effects6,  finite  density  gradients6, 
mass  ablative  effects2,  fire-polishing7,  mass  transfer8,  finite  thermal 
conduction8  and  conqpressibility8  effects. 

In  Fig.  1  we  coqpare  the  results  of  the  thin  layer  model6,  the  exponen¬ 
tial  density  gradient  model6  and  a  modified  version  (the  latent  heat  is  the 
entropy  change  across  the  ablation  layer)  of  the  boiling  film  model8  to  the 
classical  and  numerically  obtained  results.  The  boiling  film  model  includes 
mass  and  thermal  transfer  across  the  interface.  The  numerical  results  were 
obtained  with  our  two-dimensional  Cartesian  simulation  code  FAST2D4' 8  using 
a  sliding,  varibly-spaced,  Eulerian  grid.  The  sinulations  were  done  for  a  20 
in  thick,  plastic  (CH)  foil  irradiated  with  a  T.06  urn  laser  with  an  absorbed 
intensity  of  10 13  w/cm2.  The  initial  steady-state  density  profile  had  its 


density  peak  perturbed  with  a  single  sinusoidal  mode  corresponding  to  a  total 


mass  perturbation  of  1%. 
Manuscript  approved  June  2,  1983. 


The  number  of  e-foldings  was  8.0  4  0.4  for 
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all  the  modes  investigated.  It  is  apparent  that  the  classical  solution  is 

the  only  one  that  yields  the  same  k -dependence  as  the  numerical  results, 

albeit  with  a  growth  rate  larger  than  the  numerical  value. 

Recently  it  has  been  shown  that  conpressibility  is  a  destabilizing 

mechanism1 0  and  mass  ablation,  or  tip  erosion2,  is  viable  for  stabilization 

only  if  the  mass  ablation  rate  (mgb>  is  greater  than  the  rate  of  mass  flow 

into  the  spike  We  now  show  that  this  is  true  only  for  extremely  short 

wavelength  laser  light,  m  (gm  cm'^s*1)  »  pvL  *  p  v  L  where  the  subscript  c 

ar>  c  c 

denotes  the  density  (p)  and  velocity  (v)  at  the  critical  surface  and  L  is 

some  length  along  the  ablation  surface,  m  «  pv  L  where  v  is  the  velocity 

rt  s  s 

of  fall  of  the  spike.  In  the  linear  regime,  the  anplitude  of  the  spike  is 

a  =*  aQcosh  ((kg)1/2t);  thus  for  short  times,  say  one  e-folding, 

vg  *  0.3(X^g)1/2  (for  a^  *  O.IX^,  ^rt  *  2irA»  the  perturbation  wavelength). 

Now  g  *  P  /M  where  P,  is  the  ablation  pressure  ( dynes-cm"2 )  and  M  *  pa  is 

cl  a 

the  mass  (gm-cm-2)  of  the  spike  of  anplitude  a.  In  the  linear  regime  a  < 

1/2 Xrt  and  for  short  times  a  ■  0, l5Xrt.  Prom  steady  3tate  planar 
ablative  flew  models11  we  have  v_  *  9.3  x  102(IX2)1/3  cm/sec  and  Pa  *  5.6 
x  103(I/X)2'3  dynes/cm2,  where  I  is  the  absorbed  irradiance  (W/cm2)  and  X  is 
the  laser  wavelength  (urn).  This  gives  a  condition  on  laser  wavelength  for 
mass  ablation  to  dominate  RT  growth  that  X(  un)  <  3.0  x  io~2  A/ ( Z+ 1 ) p”1 /2, 
where  A  is  the  atomic  number,  Z  is  the  charge  state  and  p  is  the  density  of 
the  slab  (gm-cra"3).  Note  that  for  this  model  the  laser  intensity  and 
perturbation  wavelength  scale  out  of  the  problem.  When  the  spike  amplitude 
is  about  half  the  perturbation  wavelength  the  spike  goes  into  freefall5. 
under  those  conditions  the  spike  velocity  can  be  written  as  vs  *> 

<2Pa/p)1/2  »nd  X  <  1.6  *  10“2a/ (Z+1) P-1 /2;  thus,  a  more  stringent 
condition  applies  in  the  nonlinear  regime.  Since  the  above  models  have 


failed  to  adequately  explain  the  reduced  linear  RT  growth  and  the  k1/2 
dependence  we  propose  here  a  new  ablative-like  inhibition  mechanism  -  the 
ablation  of  vorticity. 

In  the  case  of  a  finite-sized  density  gradient,  the  discontinuous  vortex 

sheet  of  the  classical  RT  problem  is  replaced  by  a  finite-thickness  vortex 

-*■  _  + 

layer.  The  rate  of  change  of  the  vorticity  ( u  *  us  *  V  x  u(x,y))  following 

z 

♦  ♦  ♦  ♦  A  ♦  ♦ 

a  fluid  element  is  given  by  dw/dt  *  dui/9t  +  u*Vu>  «  Vp  x  Vp/p*  -  where  p 

is  the  mass  density  and  p  is  the  pressure.  The  second  term  on  the  right^iand 
side  is  the  coiqsressibility  term  and  is  typically  quite  small  in  our  simula¬ 
tions.  The  first  term  on  the  RHS  is  the  baroclinic  vector  and  is  the  source 

term  for  the  vorticity.  Fig.  2a  shows  the  contours  of  constant  du/dt 
for  a  20  urn  wavelength  1%  density  perturbation. 

The  vorticity  that  arises  from  the  differential  twisting  of  the  density 
and  pressure  surfaces  feeds  the  "heavy"  material  into  the  spike  and  the 
"licfat"  material  into  the  bubble.  Ihrt  of  the  vorticity  generated  in  the 
unstable  layer  is  ablated  and  periodically  shed  into  the  blowoff  thus 
reducing  the  linear  RT  growth.  Figure  2b  illustrates  this  shedding  when  the 
RT  instability  is  in  the  linear  regime  (3.97  ns).  The  portion  of  the  foil 
which  is  convex  towards  the  laser  (the  spike)  is  now  shedding  vorticity  into 
the  blowoff.  As  the  magnitude  of  this  vorticity  increases  it  draws  fluid  of 
opposite  vorticity  across  the  wake  cutting  off  the  growing  vortex  from  a 
further  supply  of  circulation.  The  vortex  then  separates  from  the  ablation 
layer  and  is  convected  downstream  into  the  blowoff.  The  interaction  of 
vortices  of  opposite  sign  results  in  a  relatively  constant  periodic  shedding 
from  both  the  convex  (spike)  and  concave  (bubble)  portions  of  the  foil  while 


Che  system  is  in  Che  linear  regime.  This  vortex  shedding  process  is  very 
similar  to  Che  £lov  about  a  bluff  body12  characCerized  by  Che  universal 
Strouhal  number1 3  a  dimensionless  shedding  frequency.  (The  baroclinic 


term  plays  Che  role  of  Che  bluff  body  viscous  boundary  layer).  S 


w 

where  f  is  the  shedding  frequency,  dw  is  Che  width  of  Che  wake  and  uw  is 
the  velocity  at  separation.  S  ■  0.18  and  is  constant  for  a  wide  range  of 
Reynolds  nunbers  and  body  shapes.  We  compare  the  measured  (between 
successive  sheddings  from  the  spike)  and  calculated  (f  *  0.18  uw/dw) 
shedding  frequency  as  a  function  of  perturbation  wavelength  in  Fig.  3.  To 
within  a  constant  factor,  the  agreement  is  quite  good  for  the  longer 
wavelength  inodes.  The  anomalous  increase  in  the  measured  shedding  rate  for 
the  short  wavelength  modes  is  similar  to  the  rapid  increase  in  Strouhal 
frequency  when  pairs  of  cylinders  are  close  together. 14 

To  account  for  the  reduced  linear  RT  growth  as  a  result  of  this  vortex 
shedding  process  we  write  the  RT  growth  rate  as  Yvs  *  (kg)  1/2  -  of,  where 
a  is  some  fraction  of  the  measured  shedding  frequency.  We  set  o  ■  0.4  to  fit 
the  20  pm  mode.  Figure  1  shows  good  agreement  between  the  computational 
growth  rates  and  the  rates  obtained  from  the  vortex  shedding  theory  (X).  An 
a  ■  0.4  corresponds  to  a  defect  ratio  of  0.4;  i.e.  a  measured  loss  of 
vorticity  due  to  vortex  merging  at  the  ablation  layer  of  approximately  60Z 
and  is  in  agreement  with  the  analogous  flow  past  a  bluff  body12. 

As  the  spike  grows  in  amplitude  it  begins  to  shield  the  bubble  from  the 
ablation  process.  As  a  result,  the  vorticity  is  advected  down  the  sides  of 
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Che  spikes  until  Che  temperature  gradient  is  sufficiently  large  to  cause 
ablation  (Fig.  2c).  The  linear  RT  growth  begins  to  saturate  when  the  ampli¬ 
tude  of  the  spike  is  equal  to  approximately  half  the  wavelength  of  the 
perturbation.  The  sides  of  the  spikes  cease  ablating  and  the  vorticies  be^m 
collecting  behind  the  heads  of  the  spike  causing  the  KH  roll-up  of  the  spike 
tips  proportional  to  the  wavenunber  reported  previously4’ 5  (Fig.  2d). 

It  is  well  known14  that  when  the  separation  between  pairs  of  cylinders 
is  less  than  twice  the  diameter  of  the  cylinder  the  Strouhal  frequency 
increases  well  above  that  for  a  single  cylinder  by  a  factor  of  -  5/3.  We 
find  a  similar  increase  in  Che  shedding  frequencies  for  the  short  wavelength 
cases.  This  leads  to  a  scaling  law  that  relates  the  separation  of  the 
vortices  (Aft/2)  to  their  area  and  predicts  an  increase  in  the  Strouhal 
frequency,  thus  enforcing  the  short  wavelength  cutoff  of  the  RT  mode,  for 
Art  <  8b  where  b  is  the  thickness  of  the  shear  layer  at  Che  ablation 
front,  b  scales  as  I”1 1 3  and  is  independent  of  the  foil  thickness.  Figure  4 
compares  the  predicted  short  wavelength  cutoff  value  (vertical  arrows)  to  the 
FAST2D  sioulation  results.  Also  see  Fig.  1. 

Vortex  dynamics  may  also  be  used  to  explain  the  results  of  dynamic 
stabilization  of  long  wavelength  modes  reported  several  years  ago  15.  It  is 
also  well  known  that  if  a  bluff  body  is  oscillated  in  line  with  or  transverse 
to  a  uniform  flow  with  a  frequency  near  the  Strouhal  shedding  frequency,  the 
vortices  will  shed  at  the  impressed  frequency1 3.  In  the  ablation  case  the 
equivalent  effect15  was  achieved  by  oscillating  the  laser  intensity  by  25Z  at 
a  frequency  of  3xl09s-1.  This  would  give  a  theoretical  growth  rate 
of  Y^  *  (kg)1/2  -,4fda  *  3.0xl088-1( (kg) 1/2  •  l.5*109  s”1),  considerably 


5 


reduced  below  the  unstabilized  computational  growth  rate  for  the  100  im  mode 


of  Yc  -  9.8  x  108  s"1,  in  agreement  with  the  computational  results. 

Another  important  consequence  of  the  vortex  dynamics  is  that  it  permits 
an  analogy  with  Che  magnetic  field  structure.  The  baroclinic  source  term 
also  contributes  to  the  growth  of  a  dc  magnetic  field18  and,  assuming  the 
compressibility  term  is  small,  we  find  B(Gauss)  ■  -1.04  x  10“H  w,  where 

A  is  the  atomic  number,  Z  the  charge  state  and  u>  the  vorticity  (s”1).  Using 
this  analogy  we  find  typical  field  strengths  at  the  ablation  layer  for  all 
Che  RT  modes  investigated  of  20-30  kG  in  the  linear  regime  (after  about  2  e- 
foldings)  and  3-4  MG  at  the  saturation  level  of  the  RT  instability.  The 
magnetic  field  saturates  at  the  same  time  as  the  RT  growth.  These  are 
toroidal  magnetic  fields  and  as  a  result  of  the  ablation  process  and  the  fact 
that  the  magnetic  Reynold's  number  R^  >  1  they  fill  the  overdense  region. 
Magnetic  field  generation  due  to  the  RT  instability  has  also  been  similated 
with  a  magnetohydrodynamic  particle  code17,  however,  the  profiles  and  initial 
conditions  in  that  model  were  not  characteristic  of  Che  typical  laser-abla¬ 
tive  experiments. 

Magnetic  fields  will  be  generated  at  the  ablation  layer  whenever  the 
density  and  pressure  profiles  become  non-col linear ;  i.e. ,  by  the  RT  instabil¬ 
ity,  laser  asymmetries  and/or  target  imperfections,  which  may  explain  thermal 
flux  inhibition.  For  a  series  of  2:1  laser  asymmetry  runs  for  various  inten¬ 
sities  we  find  that  B  tails  off  as  r”1  away  from  the  ablation  layer  and 
reaches  a  constant  average  value  of  <B>  ~  10-15  kG  out  to  the  critical  sur¬ 
face  with  a  resistive  decay  time  3  times  longer  than  the  3.5  ns  pulse. 


idbaifiafew-A 


By  comparing  the  thermal  flux  across  a  magnetic  field  (q^)1®  to  the  uninhib¬ 
ited  free-streaning  flux  (q^g)  it  is  easy  to  show  •  5/z(re/Xnfp)(re/L)qfs* 
where  rg  is  the  electron  gyroradius,  is  the  electron  mean  free  path  and 

L  is  the  thermal  gradient  scale  length.  From  the  asymmetry  runs  B  »  15  kG  at 
critical  and  for  1  <  6*1013  W/cm2  for  X  ■  1.06  urn  and  I  <  1011*  W/cm2  for  X  ■ 
0.53  urn,  ^  1  and  classical  thermal  flux  applies.  However, 

for  I  *  lxlO14  W/cm2  for  X  *  1.06  tin,  xj  X^f  *  0.7  and  L  *  700  un, 
thus  q^  *  .06  q^g  giving  rise  to  a  moderate  flux  inhibition  at  this 
intensity.  Assisting  the  same  B  *  15  kG  and  scaling  Tg  o  I2/5  and  L  a  I2/5, 
then  we  would  predict  a  flux  inhibition  of  0.01  for  1.06  urn  laser  light  at 
101 5  W/cm2.  This  flux  inhibition  factor  is  somewhat  smaller  than  those 
proposed  to  explain  experimental  results.  This  is  most  likely  due  to  our 
overestimate  of  the  magnetic  field  at  critical  since  we  have  dropped  the 
compressibility  term.  However,  the  flux  inhibition  from  the  magnetic  field 
structure  scales  properly  with  laser  intensity  and  wavelength. 

In  conclusion,  we  have  shown  chat  the  vortex  dynamics  of  laser  ablation 
explains  the  inhibited  RT  growth,  the  KH  rollup,  the  short-wavelength  cutoff, 
dynamic  stabilization  and  provides  a  possible  mechanism  for  thermal  flux 
inhibition  resulting  from  magnetic  fields  in  the  overdense  region. 
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Figure  1 

Comparison  of  Che  computational  linear  growth  rate  ( — * — )  to  the 

classical  solution  (vc*  - )  the  thin  layer  solution  - )  the 

exponential  gradient  solution  Yex»  - )  and  the  mass  and  thermal 

transfer  solution  (Tmt * - ).  The  X' s  are  the  theoretially 

predicted  growth  rates  and  the  arrow  indicates  the  predicted  short  wavelength 
cutoff. 
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Figure  3 

Comparison  of  Che  measured  shedding  frequency  (x)  Co  Che  value  calculaCed 
wich  a  conscanc  universal  Scrouhal  frequency  of  0.18  (0). 
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Figure  4 

Comparison  of  Che  predicted  shore  wavelength  cutoff  value  (vertical 
arrow)  to  the  simulation  results  for  7um  C,  1013  W/ca2  (0),  8  un  CH,  101  3 
W/cm2  (x)  and  7jsn  C,  1011*  W/cm2  (O). 
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